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Abstract 

It is sketched how a monostable rf- or dc-SQUID can mediate an inductive couphng between 
two adjacent flux qubits. The nontrivial dependence of the SQUID's susceptibihty on external 
flux makes it possible to continuously tune the induced coupling from antiferromagnetic (AF) 
to ferromagnetic (FM). In particular, for suitable parameters, the induced FM coupling can be 
sufficiently large to overcome any possible direct AF inductive coupling between the qubits. 

The main features follow from a classical analysis of the multi-qubit potential. A fully quantum 
treatment yields similar results, but with a modified expression for the SQUID susceptibility. 

Since the latter is exact, it can also be used to evaluate the susceptibility — or, equivalently, 
energy-level curvature — of an isolated rf-SQUID for larger shielding and at degenerate flux bias, 
i.e., a (bistable) qubit. The result is compared to the standard two-level (pseudospin) treatment 
of the anticrossing, and the ensuing conclusions are verified numerically. 
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I. INTRODUCTION 



Superconducting devices made of mesoscopic Josephson junctions are one of the few re- 
maining viable candidates for the reahzation of quantum bits (qubits). The major types 
encode quantum information in the charge, phase, or flux degree of freedom Quantum 
coherence in such devices was conclusively demonstrated at the turn of the century. While 
much remains to be done in improving the coherence times, control, and readout of individ- 
ual qubits, another key challenge on the road towards quantum-register prototypes is the 
coupling of several qubits. In the past years, coupling was achieved for each of the mentioned 
qubit types 0, 

In all published experiments, the coupling strength is a constant determined at the time 
of fabrication. For each model of quantum computing (gate-operations, adiabatic/ground- 
state, etc.), however, it would be very desirable to have the coupling tunable in situ, even if 
this requirement can sometimes be circumvented in "software" In particular, whenever 
the coupling can be continuously tuned between two values with opposite signs, it can be 
switched off naturally. 

For charge qubits, a design fulfilling these requirements was presented in Ref. j^. The 
coupling between the two qubits is mediated via a third device, the bias parameters of which 
can be used to tune the coupling constant. While the coupling element thus is nontrivial, it 
is important that it be passive, i.e., that it does not get entangled with the qubits so that 
the system's total effective Hilbert space remains four-dimensional. In practice, this means 
that the coupler must stay in its ground state, adiabatically following the qubit dynamics. 

Mediated coupling can be adapted to flux qubits 0. This paper's purpose is to contribute 
to its analytical theory. The main results can be obtained from a purely classical analysis 
of the system's potential, which may have some independent pedagogical interest. This is 
carried out in Sec. |Tllfor two choices of the coupler, both SQUID devices. A full quantum 
analysis is presented in Sec. IIIII the main result for the tunable coupling constant is highly 
analogous to its classical counterpart, but some additional effects are uncovered concerning 
a renormalization of the qubit parameters. In the formula for the coupling constant, the 
nontrivial factor is an exact expression for the coupler's ground-state magnetic susceptibility. 
In Sec. IIVI this is applied to a different regime, where the coupler behaves like a qubit itself, 
allowing for a comparison with the predictions of the two-state model. Some concluding 
remarks are made in Sec. 



II. CLASSICAL ANALYSIS 
A. rf-SQUID coupler 

Consider three rf-SQUIDs a through c in a row (Fig. Q). The a- and c-devices are 
supposed to work as (effectively) degenerately biased flux qubits, while the 6-SQUID is a 
coupling element, tunable by an external flux bias. To derive the indirect a-c coupling it 
suffices to consider the potential, consisting of both Josephson and magnetic terms {h = 1 
throughout), 

U = -EaCOs4>a - EhCOs4>b - COS 0, + (0-(/)")TL-1 (0- 0^) , (l) 
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FIG. 1: A system of three rf-SQUIDs, of which the two outer ones function as qubits and the inner 
one as a coupler. 
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That is, the direct AF inductive a-c couphng is ignored. In principle, this can always be 
achieved through a gradiometric layout. However, as long as Mac is small, one can probably 
just as well simply add an interaction Maclah to the final result Q below. In accordance 
with the roles played by the various loops, the bias fluxes (in phase units) are chosen as 

/7r + (M,,/L,)(0[°)-0^)\ 
r=[ ^b . (3) 

V7r + (M,,/L,)(0(°)-0^)/ 

Thus, (pac compensate externally for the shielding flux which the 6-loop couples into the 

a- and c-SQUIDs; see © for the precise definition of (f)^^\ 

The 6-SQUID should act as a passive coupler without introducing additional degrees of 
freedom, so bistability must be avoided. This can be achieved either by biasing it with a flux 
close to an integer number of quanta (0^ small), or for any flux bias by taking the shielding 
coefficient jSb = Ae'^LfjEi, < 1 (in SI units, f3 = 2nLIc/^o). The calculation that follows is 
valid in either case, and their relative merits will be discussed afterwards.^ 

We proceed by expansion in M. Without suggesting that the regime of M/L ^ 1 
(say, distant loops) is the most practical, this leads to a transparent result, which can guide 
numerical studies in the general case. The 6-mediated a-c coupling is expected to be O(M^), 
so the junction phases will be written as 



(4) 



^ A /3b < 1 can clearly be achieved by reducing E}, while keeping the area of the 6-loop, and hence the 
inductive coupling, appreciable. If required, a shunting capacitance could be placed across the 6-junction 
to keep the 5-SQUID in the flux regime, where it has the largest response. Also, even a SQUID with a 
unique potential minimum will have excited states, corresponding to plasma oscillations uj-p. While these 
are outside the present scope, the corresponding excitation energies should exceed any transitions in the 
a- and c-devices if the &-SQUID is to remain passive. There is an apparent conflict, since the shunting 
capacitance contemplated above would lower Wp; however, the relevant energy scale, given by the tunnel 
splittings in qubits a and c, is so small that all requirements can be met simultaneously in practice. 
Section Hill has a few more details on the system's quantum aspects. 
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these can be determined by solving the equihbrium condition V^t/ = 0. In leading order, 
the 0^°'' are just the stationary phases of an isolated rf-SQUID, 

/3,sin0f + 0f -7r = (j = a,c), (5) 
Asin0[°) + 0(°)-0^ = O, (6) 

where (j)f}c — tt can have either sign. To the same order as one has 

+ MlJLlL, Mab/LaU MMc/LaUL, \ 

= Mab/LaLt 1/Lb + MlJL^Ll + Ml/LlL^ MbJ UL, + ^(M^) . 
\ MabMbJ LaLbL, Mbc/LbL, l/L, + Ml/L,Ll J 

(7) 

While (|7j) may look tedious, using it consistently actually leads to significant cancellation 
below. 

In 0{M), one finds that 0i^c = due to the special choice of (p^ in while the 
6-loop picks up the shielding fiuxes of the neighboring ones, with a prefactor reflecting its 
susceptibility: 



I + Pb cos - 



'^(.-0i°)) + ^(vr-0f) 



(8) 



Calculating the (p^^^ in (9(M^) is straightforward but unnecessary: since one expands around 

a minimum of f/, the 0^- do not contribute to the relevant order. 

All that remains is to substitute ©^(0) into f/ in (P). Since, e.g., (01°''— vr)^ does not 
depend on the qubit state, one is left with 



Ae^LaLbLc 1 + /5feCOS(/.fe 
(3b cos (^^^^ 
L,(l + /?,cos0r) 



const + MabMbc ^ ^ (,)^ Uc . (9) 



The product of mutual inductances is expected geometrically, while lah oc cr^cr^ in qubit 
parlance. The fraction is worth a closer look: 

Pb cos ^ j pbcoscpl, /?b < 1 , ^^Q^ 



l + ^feCOS0[°^ l^' /?b>l- 

Thus, for small 0^, the coupling is AF, but it changes sign to FM as 0^ ^ vr (only attainable 
for Pb<l). 

This behaviour, and the coupling mechanism in general, have a transparent interpretation. 
Consider the self-flux curve 4>b^\(j)b)~4'b ^^e free 6-SQUID (Fig.|2)), which up to constants 
simply is the shielding current versus applied flux. If the a-qubit flips from t to |, then due 
to the mutual inductance —Mab < this effectively increases 0^. For small 0^, the curve 
has slope < 0, so that this increase generates a self-flux in the J,-direction. In its turn, the 
latter, through —Mbc < 0, acts as an |-bias for the c-loop, favouring the j-state there, i.e., 
opposite to the state of the a-qubit. This also explains the bound of unity exemplifled by 
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FIG. 2: The self-flux 0^ U'Afe) - for an rf-SQUID coupler with fib = 0.95, according to Eq. ®. 

the second line of ()10|) . since the maximum AF response is perfect shielding for (3h ^ oo (an 
uninterrupted superconducting loop). 

However, near (pl = tt the argument works in the opposite direction. There, the self-flux 
curve has slope > so that, differentially, the self-flux does not shield at all but actually 
cooperates with the external flux increase. Thus, the coupling changes sign to FM. Moreover, 
as I 1, the slope of this curve increases without bound (a precursor to bistability), so 
that the FM coupling can actually be large; this corresponds to the zero in the denominator 
of ()10|) . On the one hand, the divergence is never realized in practice, since for any finite 
M's one deals with finite differences, not slopes, on this curve (a nice counterpart to the 
discussion in for charge devices^). On the other, this potentially large FM coupling makes 
it possible to overcome any residual direct AF coupling through Mac- Since (fTUj) shows that 
the large-/3fe regime is quite inflexible, and that for -C 1 the coupling strength is limited 
by the small shielding flux, this case /?b < 1 is thought to be the preferred one. 

The above is mostly a straightforward derivation from (P), so very few specifics of the 
circuit are involved. For instance, the a- and c-SQUIDs could also be placed inside a large 
6-loop, which changes the signs of both Mab and M^c, so that Q is invariant. This design 
modification makes it clearer that the 6- loop is mostly a fiux transformer 0], with a weak 
link providing the tunability. Moreover, the final result (jH} depends only trivially on the 
properties of the a- and c-devices — through the fiux which these couple into the 6-SQUID. 
Hence, the generalization to other types of fiux qubits should be obvious. 

B. dc-SQUID coupler 

The above indicates that the simple rf-SQUID should be a satisfactory coupling element. 
Still, it is worthwhile to also examine a dc-SQUID j^], which provides additional fiexibility 



^ Indeed, since the shielding current is the derivative of the energy with respect to the external flux, the 
induced coupling strength is proportional to the second derivative or finite difference of the SQUID band 
structure, in complete analogy with j^. 
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FIG. 3: A dc-SQUID mediating a coupling between two rf-SQUID qubits. 



since it can be both flux- and current-biased. In particular, if only a current bias would 
turn out to suffice for coupling tunability, that could be preferable in situations where a ffux 
bias is problematic due to stray fields etc. The 6-device now has separate left and right arms, 
designated as 61 and 62, together carrying (Fig. El). Again ignoring the direct a-c mutual 
inductance (not that this makes much difference anywhere), we thus take the coil-to-coil 
inductance matrix as 

/ L„ -M„,i -Mab2 \ 

Lb2 
-Mb2c 



\ 



-Mabl 
-Mab2 





-Mabl 

Lbi 
Mbu 
-Mbic 





-Mbic 

-Mb2c 
Lr 



(11) 



However, always — hi + ffux-current equation relates three- vectors pertaining 

to the loops: 

$ - = L/ , (12) 

with 



ab2, Lb2, —Mb2c) ^b 

T 



$ = ($a,$b,$ 

= + (-M, 
I 

La -Mab 

L = I -Mab Lb -Mbc 
-Mbc Lc 



M, 



ab 



M, 



abl 



M, 



ab2 



Mbc = Mblc + Mb2c , 



Lbi + L 



b2 



^bl 



Lbi + Ml 



bl2 



^b2 



Lb2 + Ml 



612 



(13) 
(14) 
(15) 

(16) 

(17) 
(18) 



On the other hand, the four junction phases are all independent dynamical variables. They 



are related to the ffuxes by 



%l+4>b2, 



2e$. 



One should first of all determine the Hamiltonian. Of the Kirchhoff (circuit) equations, 
the first four are simply the Josephson relations 0^ = 2eQa/Ca etc., with Qa the junction 
charge and Ca its capacitance. The other four express current conservation: 



Qa 
Qbl 

Q 



b2 



-2eEaSm4>a - h = -2eE,sin0„ - [L-i($-^")]„ , 
-2eEbi sin06i-[L-i($-$^)], , 

-2eEb2 sin 062 - 41 - lb = -2e^fe2 sin 0^2 - [L-i($-$^)], - J? , 



Qc = -2eEcSm(f)c - [L-^($-$^) 



(19) 
(20) 
(21) 
(22) 
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It is readily seen that these are the Hamilton equations 
-2edH/d^j, for 



2edH/dQj and Qj 



H 



E 



Qi 

2C, 



Ej cos (f)j 



2e 



(23) 



with j G {a, 61, 62, c}. It is less apparent that H is actually symmetric in the two SQUID 
arms, as it should be: if is a function of the junction phases and charges only, and the 
asymmetric choice of Ibi as the loop current was merely an interim convention. That is, 
designating the last two terms of ()23|) as "magnetic" and "bias" energy, respectively, is a bit 
arbitrary. Working out L~^, one sees that up to a constant, one can rewrite 



H 



j 



(24) 



where 



2e det L 



i^bias = [{Mab2Lbl-MablLb2)L,+ {MablMb2c-Mab2Mblc)Mbc](j)a 
+ [Mab2MabL, + LaMb2cMbc - L^Lfes^eJ^W 
+ [-MablMabL, - LaMblcMbc + L,LwL,](/.62 
+ [La{LblMb2c-Lb2Mblc) + Mab{Mab2Mblc-MablMb2c)]<Pc 



(25) 



with det L = LaLbL^ — LaM^^ — M^^Lc. It would be instructive to see if the above, especially 
the form for iibias, can be reproduced (presumably almost mechanically) by the "node 
formalism" of circuit analysis l^] ■ 

From here on, the analysis is a straightforward generalization of the rf-SQUID case, again 
writing the phases as in (jU, and expanding iibias in to the same order. We again want 
the a- and c-qubits to be effectively at degeneracy, which presently means that their own 
external flux should compensate also for J^: 



TT 



0^ 



L-b'[MMfh€-<K] 



2eI^{Mab2Lbl-MablLb2)T 



(26) 



,7r + L,-i[M,,(0i?+0g^-0^) + 2eIl{LbiMb2c-Lb2Mbic)]. 
In (9(M°), the stationary phases obey the standard equations for the isolated devices: 

' (j = a,c), (27) 

, (28) 

(29) 



f3bi sin (t)f^ + 



(3b2 sin (pf^ + 



(3j sin0^"' + (py - TT 
(0) _ 2eRU 



. 



Here, [3a{c) = ^^^La{c)Ea{c) is standard, but note that the definitions [3bj = ie^LbEbj (j = 1, 2) 
involve the full loop inductance of the 6-SQUID, not the individual arm inductances. 

In 0{M), our choice ()26|) for 0^ again ensures that (f)^a^ = (j)^c^ = 0, while the equations 



for 



4l(2) 



are coupled, as were (PHj) - (f^ for 



course linear): 



.(0) 
\l{2) 



(although, unlike the latter, they are of 




l+(3bi cos 0^°^ 



l+f3b2 cos 0^2^ 



Mab, (o)n , Mbc, 



'>(0)^ 



(30) 



Like for the rf-SQUID, it is best to leave the 0^- unevaluated in the M-expansion of U, 
and one sees that their contribution cancels (to this order), since one expands around a 
potential minimum. Taking advantage of some further cancellation and dropping all terms 
which do not depend on the qubit state, one arrives at 

TT f , /1# /3fei/362COS0['J^COS0[J 

U = const + MMc — joy— ,(0)^^ ^ -70^ T(5k ^-^'^ ' (^l) 

LbiPbl cos 0^/ + fJb2 COS + PblPb2 COS (pl^' COS ) 

in which all tunability thus comes via the dependence of ^^^jg) on (06, ^b) given in 
fl28|l - (j^n|l . For a simple consistency check, let us take — » and Eb2 oo. Then, 
fl^ shows that 0^2'' = 0; while (^Hjl reduces to jSbisiiKpf^ + (f)f^ — 0^ = 0, and the large 
fraction in becomes (3bi cos^^'J''/ [l + jSbi cos^^'J''] ; in view of the remark below ([2^1) . this 
is exactly the result for rf-SQUID coupling. 

Another well-known case is the symmetric dc-SQUID with negligible shielding, for which 
(EHl)-® reduce to 0[f + 0[J = 0^ and 2/c cos(0^/2) sin(0[°^-0^/2) = I^. Even in this 
simple case one can have FM coupling, for cos^(0^/2) < |/^|/2/c < |cos(0^/2)|. While one 
thus needs a nonzero flux bias, tunability of suffices for changing the sign of the coupling. 

One notes a few features, and it should be investigated whether these carry over to more 
generic devices described by our equations: (a) to achieve FM coupling, one needs nonzero 
flux and current bias; (b) the denominator in is always positive; (c) cos^^'J'' and cos^^j'* 
never become negative simultaneously. Need for current bias: for a symmetric SQUID this 
is easy to prove, since for = one then has ipf^ = (j)^^2 y ^^Y 0b there is a stationary 

state with |0['j'''| < 7r/2- This cannot be generalized further, since we have already seen that 
the (asymmetric) dc-SQUID contains the rf-SQUID as a special case, and the latter can 
mediate FM coupling without any current bias. Need for flux bias: again easily proved for 
symmetric devices with arbitrary shielding, which have 0^°^ = — 0^2'' etc. for 0^ = 0; in 
general, however, ((211)^(1211) show that an inductance imbalance between the two SQUID 
arms can play the same role as a nonzero 0^. Finally, note that the rhs in (jHUjl features the 
second-derivative matrix for the potential of the free 6-SQUID (characterizing the device's 
susceptibility). For a stable minimum, this matrix should be positive-definite, from which 
(b) and (c) immediately follow in general. (As remarked above, the second derivatives keep 
appearing because of the M-expansion. For finite M, one has second differences instead j^.) 

One of the possible uses for a dc-SQUID coupler could be to moonlight as a readout device 
through switching-current measurement, so let us treat that in the present notation. Rather 
than think of the maximum bias for which (j28p - (j29p have solutions, it is convenient to fix 0^ 
and consider as a function of the junction phases, which are still related among themselves 
by one nontrivial constraint. Critical bias now corresponds to dl^/dcf) = 0; working things 
out readily yields 

f3bl cos 0^°^ + (3b2 cos 0^2'' + (3blf3b2 COS 0^°^ COS 0^2^ = (32) 

as the relevant condition.^ To calculate switching- current sensitivity, we now vary 0^, and 



^ Compare and the end of the previous paragraph: one can equivalently characterize critical bias as the 
current for which the potential becomes marginally unstable. We note in passing that, for near-critical 
bias, the induced coupling is always FM, except for the doubly nongeneric points given by ()34|l and H35|) . 
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see which variation in conserves the above relation. One finds 



7 TX,max n i(0) 

_ /36lCOS0^/ _ ^^^^ 



#6 2e[Lfe + L6i/3wcos0[?] 



using (|32p . one again finds that this is more symmetric in the SQUID arms than it looks. 
So the sensitivity vanishes iff, at critical bias, both junctions become critical individually. 
This happens for 

" ^ ' " ^ 2eL, • ^^^^ 

As long as the "critical fiux" product Lbjlc,bj is different for the two arms, one has a finite 
sensitivity at 0^ = and in fact for any 0^, since one can ramp with either sign. For the 
sake of completeness we also give the zero-sensitivity points 



± 



LbiPbi + Lb2f3b2 

TT 

Lb 



,Pb2-Pbl /oc\ 



which correspond to minima in the switching current. 

Finally, the above analysis only shows that the potential part of ()24|l - (j^K|l enables tunable 
coupling in the classical case. That is, there could be devices outside its regime of validity 
which are nonetheless suitable couplers, as long as the transformer remains passive, i.e., 
confined to its lowest energy state/band j^. For a quantum analysis, the interacting, biased 
Hamiltonian (j^^ - ()25|) can be used as a starting point. The method of M-expansion should 
again reduce the problem to one for the uncoupled, biased, dc-SQUID. Below, however, the 
quantum case is taken up for rf-SQUIDs only. 



III. QUANTUM ANALYSIS 

Given that quantum tunneling is crucial for the qubits a and c, a consistent quantum 
treatment of the whole system is desirable. In the quantum case, the M-expansion works 
as above. The other crucial step is an adiabatic approximation for the 6-coupler (cf. |5|), 
replacing operators with their ground-state expectation. This very language is only mean- 
ingful if the rf-SQUID 6 is a well-defined separate entity, i.e., for weak coupling. Specifically, 
the coupling energies must be much smaller than the first excitation energy of the 6-device, 
but this does not restrict their size compared to the level splittings of qubits a and c. 

Thus, let us collect all terms in the Hamiltonian H = J2j=ab cQ^'^^j ~ ^ [with U as 
in ([Q)] involving the coupler. To the relevant order in M one has 

^,.5L_^,eos0. + ^|ri + i^+^^^- 



2Cb -^-^^"8^u^z:ii"iii:r^^ 

+ 2^(0,-0^)(0,-0^) + 2^(0,-0^)(0,-0^)} (36) 



2Cb-^''''''^'^ 8^ 



- Eb cos (j)b + " oJreS ' (37) 
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with 



Lf = Ljl--^--^] , (38) 

<pr' = 06 + s<pt = 0^ + + . (39) 

The adiabatic step amounts to the replacement 

ft-ffW-'j-if (40) 
and for small coupling, one can further write 

r(0) (0) 

£f\<Pr^) ^ Si'\4>t) - -j^6<Pt + ■ (41) 

r(0) 

as 



Here, is the ground-state expectation of the loop current; the susceptibility is defined 



xi°'--'£-U^-'-^], (42) 



which gives the bound xi'''* < ^/Lb (no absolute value on the Ihs) — the maximum shielding 
is perfect diamagnetism for an uninterrupted loop [cf. below fjlOj) ]. 

All that remains is to substitute ()40p and ()41|) back into the full H; after some cancellation, 
one finds 

^ = 1^ - cos 0, + ^^^^J^^^ ^ + (a - c) + H,^, , (43) 

- 2e/f ^M^fe , (44) 
Ml 



Lf = lAi- "^^xT ] , (45) 



-'a 
.(0) 



H,^, = MMcXrialc . (46) 

This does not involve assumptions on the states of the a- and c-qubits, or on their flux 
biases. Note that ()46|) involves current operators for the a- and c-devices in contrast 
to, e.g., (jl^ and (jH)). Equations (jl^ and (jl^ imply that the qubits' effective parameters 
depend subtly on the coupler's working point, presenting a possible experimental challenge. 

As in the classical case, the perturbative step has reduced the problem to that of the un- 
coupled 6-SQUID, and we henceforth omit its index as well as the zero superscript indicating 
free devices. Evaluating ()42j) in the potential minimum, one verifies that reduces to 
(jH)) in the classical limit. In the quantum case, numerical differentiation in ()42j) presents no 
problems, but it will be instructive, and useful for Sec. lIV| to express x in terms of quantities 
at a single bias point. To this end, is evaluated in terms of ^o(0) = dfjj^i'o^cj)) , where 
normalization implies that (ipol^o) = 0. In its turn, is solved from the ^^-differentiated 
Schrodinger equation 



{H-So)^o+{j^-Sl,]^o = 0, (47) 
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in which i/jq figures as an inhomogeneous term, and where Sq = d(pxSo = (0'^— 0)Q/4e^L. 
Using a variation-of-constant solution to ()47|). one then derives'^ 



C 



4eU ^oi^) 



C^0l(01-(0l)o)^o(01 



(48) 



This aheady expresses considerable simplification, in that two of the integrals in a 
triple integral were written as the square of k{(f)) = J^^d4>i {4>i~{4>i)o)'ipo{4>i) = 
— d(f)i (01 — (0i)o) V'ol^i)- However, /c(0) exponentially for — oo (0 ^ — oo) 
will not be realized numerically for the former (latter) representation of k. Hence, we have 
to use both. Putting /i(C) = jff^ # F(0)/^g(0) and /2(C) = d(f)P{(j))/ij^{(j)), one has 
X = L~^[l — (C/4e'^L){Zi(0) + /2(0)}]. Now, x can be found by a single integration of the 
ODE systems 

'I' ^ ^2 ^^2 n/ ^ _^2 ^^2 

A:'(0) = (0-(0g^2(0) , ifc'(0) = (0-(0)o)V^o'(0) • (49) 
J^{-oo) = k{-oo) = [/2(oo) = A;(oo) = 

A different derivation uses that, for stationary states, the shielding current equals the 
Josephson current (readily proven formally): 

It is correct to evaluate (j^UI) with the aid of (jTTj) . but this yields unwieldy asymmetric 
expressions. Instead, we observe that a change in flux bias is a uniform translation if there 
is no Josephson potential, i.e., we set ^0 = ''70 — i'o with ('0o|'7o) = and 

(H - £o)r]o + {E sin - 4)^-0 = , (51) 

using the matching expression Sq = -E'(sin0)Q. Combination leads to 

X = 4e'E{cos<P),-4CE' j -|l^|y%0i(sin0i-(sin0i)o)^o2(0i)| . (52) 



The translation of ()52|) to ODE form is wholly analogous to ()49|) and will not be repeated. 

For small (3, (j^^ is clearly preferable to (pSjl . since the double integral yields only a small 
correction to the flrst term. Compare (j50|l . readily yielding an estimate of x for weakly 
shielded SQUIDs, to where in this case one has to evaluate the difference of two nearly 
equal terms. Conversely, for large (3, (0)q ~ 27m is approximately constant over most of 
the bias range, giving a clear advantage to ()42|1 and PH|) . which capture this behaviour in 
their flrst terms. The equivalence of the four expressions for x has been thoroughly tested 
numerically. 

Results are shown in Fig. 0] It is seen that the flnite spread of the ground-state wave 
function smooths out the classical susceptibility, in particular near 0^ = tt where the FM 
response is maximal. The flnite quantum tunneling amplitude also means that there is no 
longer a sharp transition at /3 = 1, although the small tunnel splitting for larger /3 is likely 
to violate the adiabaticity requirement. 



All wave functions are taken real w.l.o.g. 
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*>[m4.J t'lmoj 1° 1° g 1° 1° 

FIG. 4: The reduced SQUID susceptibility Lx versus flux bias <I>^ for (5 = 0.95. (a) Classical 
response, (b) Close-up of the classical (solid line) and quantum response (dashed line, for g = 5-10'^) 
for biases near ^$oi the two curves are indistinguishable elsewhere, (c) The quantum susceptibility 
at = versus capacitance, showing the approach to the classical limit. 



IV. APPLICATION: SQUID SUSCEPTIBILITY NEAR AN ANTICROSSING 
A. Analytics 

While the above was derived with monostable coupling devices in mind, it should be 
general, hence equally applicable to rf-SQUID qubits close to the anticrossing, where it can 
be compared to the predictions of the two-state model 

8f{5^'') = const - ^A(5<l>-)2 + e((5$-)2 , (53) 

with 5$^ = $^ — |$o- The bias e is defined as half the difference between the "local well 
energies", which for definiteness are taken as the ground-state energies with an artificial 
nodal condition imposed at the potential maximum. If /3 — 1 is taken fixed, then the two- 
state model should be valid asymptotically in g = 2CE/e^. Note that it incorporates 
several quantum effects short of tunneling: at finite bias the curvatures of the two wells 
change in opposite directions, affecting the zero-point energies (note 15 in Ref. jll|), and the 
anharmonicity of the wells is also accounted for. The "tunneling amplitude" A is primarily 
defined through (j^Hj) . without reliance on WKB estimates. Bearing in mind that A (5$^) 
[e(5$^)] is even [odd], one finds x*(0) = ^o^(O)" = -A"(0) - J^/A for the two-state 
approximation to with A = A(0) and Jp = — e'(0). 

These predictions make it clear that x(0) is exponentially large in g, so that in the 
evaluation of (pSjl and ()52|) . their first terms can be safely ignored. Further, (0)q = tt and 
(sin0)g = by symmetry for (p^ = vr, which will be assumed henceforth. In either formula, 
the increase in x can only come from the factor ^q'^, which is large for |0| ^ oo and near 
the barrier. However, in both formulas and both for oo and (p — oo, this is overcome 
by {■ ■ ■ }^ decaying as i/jq. Therefore, the only relevant contribution comes from the barrier 
region, where to within exponentially small corrections the factors {■ ■ ■ }^ are constant, and 
readily related to the loop current. Thus, in the regime of interest, both (j48p and (j52|) 
evaluate to 

X(0) - -?| . (54) 



A ^4'U<I>) 

where means an integral over the barrier region, which for definiteness can be taken as 
everything between the two well minima (so as not to rely on "classical turning points"). To 
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relate the integral in to A, define f{(j),£) as the solution to the Schrodinger equation 
at (continuously tunable) energy S satisfying the left boundary condition. Of course, this 
function will diverge for — > oo unless S is an eigenenergy. Its energy derivative satisfies 
a Schrodinger equation in which / figures as an inhomogeneous term [cf. (j47p for the bias 
derivative], which we solve by the variation-of-constant Ansatz d^f = fh. Assuming that 
the £^-dependence of / can be linearized for Sq < £ < Si, the significance of these functions 
is that dsfin, S) ^ [0 - iIjo{'k)]/2A h{n, So) ^ -1/2A. Solving for h, one finds 

this defines h only up to a constant, since / was defined only up to a, possibly £^-dependent, 
normalization. We now fix J^^dcp f^{(f)) = |, so that h'{(f)) ^ —C/Ae'^f'^^cj)) throughout the 
barrier region. Further, with ^min the left potential minimum, h{(l)ram,S) is of order one, 
hence negligible, and one finds /i(7r, £^o) ~ — (C/4e^) d(f) i^q'^ {(p) = — (C/8e^) J^d(f) i^q'^ {(p) . 
Combination shows that = (C/4e^) J^dcp i/jq'^ {(p) , and substitution in yields 

X(0) ^ -| . (56) 

Comparing with what was found below for x'^^i^) the two-state approximation, it 
leads one to conclude that actually A"(0) = 0. More precisely, comparing the orders of 
magnitude, it indicates that A does not have relative variations of order one in the interval 
|e| ~ A, which is all one cares about. This flatness of A in the anticrossing region is often 
assumed without much discussion, let alone proof; even the comparatively detailed WKB 
treatment in Ref. 2| is complete only for highly excited states in the local wells. Yet, the 



standard WKB expression for the tunneling amplitude is very sensitive to the barrier action 
integral, suggesting that there could be an equally strong sensitivity to bias variations. 

All in all, the above paints a positive picture of the standard two-level approximation. 
Its regime of validity should follow from |e((5$^)| <^ S2 — Si, where the rhs is of the order of 
the plasma frequency. Hence, the range of allowed 5$^ is actually algebraically decreasing 
in g so that asymptotically, variations in loop current etc. become negligible within this 
range. At the same time, the range of permissible e/A increases exponentially, so that wide 
sweeps in e, as e.g. studied in the context of Landau-Zener transitions [l^, can be performed 
without leaving this range. 



B. Numerics 

If one numerically finds the two lowest eigenf unctions of a near-degenerate bistable 
rf-SQUID, one can uniquely determine the parameters of the corresponding two-level low- 
energy Hamiltonian^ 

H^^ = -ea' - Aa" . (57) 



^ We have used similar wordings to describe the procedure in Refs. but the concepts are different. 

There, the impedance response as a function of flux bias was used to fit a single set of qubit-Hamiltonian 
coefficients. Here, we have numerical access to the wave functions, and can therefore determine these 
coefficients independently at each point in parameter space. 
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1.38 • 10-6 
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5 • 10-6 


1.38 • 10-5 


5.64450 • 10-6 
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5 • 10-5 


1.38 • 10-4 


5.65360 • 10-6 
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5 • 10-4 


1.38 • 10-3 


6.50793 • 10-6 



TABLE I: Numerically extracted biases e and tunneling amplitudes A according to H58|) . for a 
near-degenerate rf-SQUID with shielding /?, reduced capacitance g = 40, and flux bias 5^^. 



Namely, from the numerics one extracts the gap Sio = Si — Sq, plus the well distribution of 
the ground-state wave function, = J^^d(j)iljQ{(f)) and = J^d(f)'i/jQ{(j)), where the barrier 
point 0b denotes the location of the potential maximum. The average energy is an irrelevant 
constant [assumed zero in the traceless (jFfj) ]. p and q are constrained by normalization, and 
the well distribution of the first excited state is fixed by orthogonality, so that one has 
two pieces of data. Comparing these to the predictions of (|S7j), one can solve for the two 
unknowns as 

e = £io— — , A = Siopq . (58) 

Thus, the coefficients of the two-state model are extracted without any handwaving, even 
without the artifice of an extra boundary condition stipulated below (IHHjl : the only ingredient 
is a physically transparent coarse-graining of the wave functions. Note how e, A attain the 
expected limits at degeneracy p = q, and e also for large bias pq 0. 

Another way to arrive at this result is easier to generalize to multiple qubits. It involves 
using the transition matrix6 U = {p~q) (i-e., having the eigenstates for its columns) to 
rotate the Hamiltonian Hs = diag(£^o, ^i) from the energy basis back to the fiux basis, viz., 

UHeU^ = ^1 - Sj^a^ - S.opqa^ . (59) 



^ Sign convention: it is logical to have < correspond to e < 0. The former favours smaller cj), i.e., the 
left well; by H57|l . the latter favours — —1, so that the bottom row in U corresponds to the left well. 
Further, the columns in U and those in H£ of course have to be ordered consistently. 
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Results are summarized in Table HI With increasing (3, one sees the same trends predicted 
at the end of Sec. IIV Al for increasing g: the range of validity of (jFfj) with constant A 
increases, in terms of not but e/A. Remaining deviations A((5<l>^) — A(0) ~ (5$^)^ are 
too small to contradict the analytical conclusions. The breakdown of (j57p occurs when the 
wave function in the depopulated well is so small that contributions near 0b can no longer 
be neglected, precluding a coarse-grained description. Even then, the linearity of e(5$^) is 
still excellent, as shown by the last entries in Table |1] for each /3, and the breakdown is of 
little practical consequence as A is too small compared to |e| to have an appreciable effect 
on the spectrum of (jFfj) . 

The above situation for a single qubit is uncharacteristically fortunate, in that the re- 
quirement of unitarity on U has circumvented coarse-graining the excited state ipi. Already 
for two qubits, there is no natural counterpart to this shortcut, and the 4x4 matrix U 
having coarse-grained eigenstates for its columns is not exactly unitary, introducing a small 
uncertainty into the above procedure. The root cause of these subtleties is the problem- 
atic nature of the ubiquitous term "flux basis," which strictly speaking should consist of 
(5-functions in the coordinate representation. Pending further research, a completely satis- 
factory resolution may be as follows: calculate the matrix elements of in the low- energy 
subspace. This Hermitian matrix has an eigenbasis, and the eigenvectors can be taken as 
columns of a transition matrix V J Recognizing that V consists of 0-eigenvectors in the 
energy basis whereas U consists of ^^-eigenvectors in the flux basis, one should thus consider 
ifqb = V'^ H^V as an alternative to 

V. CONCLUSION 

A different sign-tunable flux transformer, relying on the gradiometric nature of the em- 
ployed qubit, was previously considered in [1^. The transformer itself is also gradiometric, 
and tunability is achieved by incorporating compound junctions with variable couplings in 
either arm. In another type of gradiometric flux transformer, each arm of the gradiometer 
couples to one adjacent device, and the tunable element is a single compound junction in 
the central leg [15|. This does not give sign-tunability, though: compare (4) in lis] to (HUD 
above. For the present type of mediated coupling, the prediction of sign-tunability is ex- 
tremely generic: for any SQUID-type device, the /'^°^($^) curve is smooth and periodic, and 
therefore will have pieces with both positive and negative slope in any $o-iiiterval. Compar- 
ing the two implementations, the rf-SQUlD (Sec. Ill A|) does not have any galvanic coupling 
to external circuitry, which should limit decoherence. On the other hand, it was already 
mentioned that a dc-SQUID coupler f Sec. Ill B]) could be useful for readout. Either seems 
promising, and we look forward to experimental results on tunable coupling, which may well 
be the next big step forward in superconducting qubit technology. 

The analysis in Sec. HVlis a prelude to the case of multiple interacting qubits, where the 
precise form of the low-energy Hamiltonian ifqb is still open to discussion. Namely, there 
could be multi-qubit tunneling, say, along the diagonals in a two-qubit potential, which 
would lead to a^-a^ and a'^-a^ interactions note in proof). There is little doubt 



^ For consistency, the columns of V should be ordered according to decreasing eigenvalues. For, say, two 
qubits, these columns will be the simultaneous eigenvectors of the commuting operators 0a, 06 in the 
low-energy subspace. 
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that such terms would be small, but this does not automatically imply that their effect 
is negligible. Namely, they could, say, cause a gap at an |||) |||) anticrossing (for AF 
coupled qubits) in first order, while conventional tunneling terms ~ contribute to this gap 
only in second order. Further, there seems to have been little investigation of the possibility 
that (single-qubit) tunneling in qubit a might depend on the state of qubit h and vice versa, 
which would lead to cr^cx^ etc. terms. To the extent that, for a-tunneling, the 6-qubit can be 
regarded as an external flux bias, the flatness of A observed in Sec. IIVI indicates that such 
terms should be very small. Finally, for more than two qubits, one can anticipate three-qubit 
etc. interactions — particularly if, for stronger coupling, one qubit 's persistent current can 
depend on the state of the others. In any case, rotation methods such as in Eq. fl^ should 
enable one to settle all these questions conclusively. 
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